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1 Introduction 

In solid-state physics, the phenomenon of localization is usually perceived as 
arising from extrinsic disorder that breaks the discrete translational invariance 
of the perfect crystal lattice. Familiar examples include the localized vibra- 
tional modes around impurities or defects in crystals and Anderson localization 
of waves in disordered media pQ. The usual perception among solid-state re- 
searchers is that, in perfect lattices excitations must be extended objects as well, 
essentially plane wave like. Such firmly entrenched perceptions were severely 
jolted since the discovery of discrete breathers (DB), also known as intrinsic 
localized modes (ILM). These states are typical excitations in perfectly peri- 
odic but strongly nonlinear systems, and are characterized by being spatially 
localized, at variance to plane wave states [3 EH]. 

DB-like excitations, being generic objects, have been observed in a large 
variety of lattice systems that include bond excitations in molecules, lattice 
vibrations and spin excitations in solids, charge flow in coupled Josephson junc- 
tions, light propagation in interacting optical waveguides, cantilever vibrations 
in michromechanical arrays, cold atom dynamics in Bose-Einstein condensates 
loaded on optical lattices, among others. They have been extensively studied, 
and a high level of understanding about their properties has been reached. 

Two decades of intensive research have polished our theoretical understand- 
ing of DBs in classical nonlinear lattices. Less is known about their quantum 
counterparts - quantum breathers (QB). This chapter is devoted to a review of 
the more recent studies in this field. The concept of QBs is closely related with 
the theme of dynamical tunneling in phase space. 
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1.1 A few facts about classical discrete breathers 



Let us study the combined effect of nonlinearity and discreteness on the spatial 
localization of a discrete breather on a basic level. For that we look into the 
dynamics of a one-dimensional chain of interacting (scalar) oscillators with the 
Hamiltonian 

H = J2 \vl + V{x n ) + W{x n -x n - X ) . (1) 

The integer n marks the lattice site number of a possibly infinite chain, and x n 
and p n are the canonically conjugated coordinate and momentum of a degree 
of freedom associated with site number n. The on-site potential V and the 
interaction potential W satisfy V'(0) = W'(Q) = 0, V"(0),W"(0) > 0. This 
choice ensures that the classical ground state x n = p n = is a minimum of the 
energy H . The equations of motion read 

X n = Pn , Pn = -V'(x n ) - W'(x n - X n -l) + W' (x n+1 - X n ) . (2) 

Let us linearize the equations of motion around the classical ground state. 
We obtain a set of linear coupled differential equations with solutions being 
small amplitude plane waves: 

x n (t) ~ e i( ^*- 9n) , co 2 q = V"(0) + 4W"(0) sin 2 (|) . (3) 

These waves are characterized by a wave number q and a corresponding fre- 
quency uj q . All allowed plane wave frequencies fill a part of the real axis which is 
coined linear spectrum. Due to the underlying lattice the frequency u> q depends 
periodically on q and its absolute value has always a finite upper bound. The 
maximum (Debye) frequency of small amplitude waves w w 

Depending on the choice of the potential V(x), tu q can be either acoustic- or 
optic-like, 1^(0) = and 1^(0) 7^ 0, respectively. In the first case the linear spec- 
trum covers the interval — < u q < which includes uj q= o = 0. In the latter 
case an additional (finite) gap opens for \u q \ below the value w = y/V"(0). 

For large amplitude excitations the linearization of the equations of motion is 
not correct anymore. Similar to the case of a single anharmonic oscillator, the 
frequency of possible time-periodic excitations will depend on the amplitude 
of the excitation, and thus may be located outside the linear spectrum. Let 
us assume that a time-periodic and spatially localized state, i.e. a discrete 
breather, x n (t + T5) = x n (t) exists as an exact solution of Eqs.(|2]) with the 
period XJ, = 27r/f2fe. Due to its time periodicity, we can expand x n (t) into a 
Fourier series 

x n (t) = Y, A ^ elkSlbt ■ ( 4 ) 

k 

The Fourier coefficients are by assumption also localized in space 

^.fc.H-j-oo -> ■ (5) 
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Inserting this ansatz into the equations of motion ([2]) and linearizing the re- 
sulting algebraic equations for Fourier coefficients in the spatial breather tails 
(where the amplitudes are by assumption small) we arrive at the following linear 
algebraic equations: 

k 2 n 2 b A kn = V"(0)A kn + W"(0)(2A kn - A fe ,„_ 1 - A k>n+1 ) . (6) 

If fcf2(, = uj q , the solution to © is A k ^ n = c\e lqn + C2^ lqn ■ Any nonzero 
(whatever small) amplitude A k ^ n wil thus oscillate without further spatial decay, 
contradicting the initial assumption. If however 

kfl b ^ LO q (7) 

for any integer k and any q, then the general solution to ^ is given by A k n = 
c\K n + C2K~ n where k is a real number depending on uj q , n b and k. It always 
admits an (actually exponential) spatial decay by choosing either c\ or C2 to be 
nonzero. In order to fulfill (J7]) for at least one real value of n, b and any integer k, 
we have to request \uj q \ to be bounded from above. That is precisely the reason 
why the spatial lattice is needed. In contrast most spatially continuous field 
equations will have linear spectra which are unbounded. That makes resonances 
of higher order harmonics of a localized excitation with the linear spectrum 
unavoidable. The nonresonance condition (JT)) is thus an (almost) necessary 
condition for obtaining a time-periodic localized state on a Hamiltonian lattice 
0- 

The performed analysis can be extended to more general classes of discrete 
lattices, including e.g. long-range interactions between sites, more degrees of 
freedom per each site, higher-dimensional lattices etc. But the resulting non- 
resonance condition (JT)) keeps its generality, illustrating the key role of discrete- 
ness and nonlinearity for the existence of discrete breathers. 

Let us show discrete breather solutions for various lattices. We start with a 
chain (fTJ) with the functions 

V(x) = x 2 +x 3 + ix 4 , W(x) = 0.1a; 2 . (8) 

The spectrum ui q is optic-like and shown in FigJTJ Discrete breather solutions 
can have frequencies Of, which are located both below and above the linear spec- 
trum. The time-reversal symmetry of © allows to search for DB displacements 
x n (t = 0) when all velocities x n (t = 0) = 0. These initial displacements are 
computed with high accuracy (see following sections) and plotted in the insets in 
FigQ] [3\ . We show solutions to two DB frequencies located above and below uj q 
- their actual values are marked with the green arrows. To each DB frequency 
we show two different spatial DB patterns - among an infinite number of other 
possibilities, as we will see below. The high-frequency DBs (f2& « 1.66) occur 
for large-amplitude, high-energy motion with adjacent particles moving out of 
phase. Low- frequency DBs (f2b ~ 1.26) occur for small-amplitude motion with 
adjacent particles moving in phase. 
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Figure 1 : The frequency versus wavenumber dependence of the linear spectrum 
for a one-dimensional chain of anharmonic oscillators with potentials ([8|). Cho- 
sen DB frequencies are marked with green arrows and lie outside the linear 
spectrum w q . Red circles indicate the oscillator displacements for a given DB 
solution, with all velocities equal to zero. Lines connecting circles are guides for 
the eye. From [3]. 



In Figl2]we show two DB solutions for a Fermi-Pasta-Ulam chain of particles 
coupled via anharmonic springs V(x) = 0,W(x) = \x 2 + \x 4 (c.f. {TJ) which 
has an acoustic-type spectrum [5]. The DB frequency is in both cases Oj, = 4.5. 
Again the displacements shown for an initial time when all velocities 

vanish. In the inset we plot the strain u n = x n — x n ~\ on a log- normal scale. 
The DB solutions are exponentially localized in space. 

Finally we show DB solutions for a two-dimensional square lattice of an- 
harmonic oscillators with nearest neighbour coupling. The equations of motion 
read 

Xi.j — k(xi-\-\.j Xi—i^j '2xi_j^ -\- k{x%^j-\-\ X%^j—\ '2xi_j^ j (^) 

with oscillator potentials V(x) = \x 2 + \x A . In FigO we plot the oscillator 
displacements with all velocities equal to zero for three different DB frequencies 
and k = 0.05 [6]. For all cases adjacent oscillators move out of phase. 

We conclude this section with emphasizing that DB solutions can be typically 
localized on a few lattice sites, regardless of the lattice dimension. Thus little 
overall coherence is needed to excite a state nearby - just a few sites have to 
oscillate coherently, the rest of the lattice does not participate strongly in the 
excitation. 
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Figure 2: Discrete breather solutions for a Fermi-Pasta- Ulam chain (see text). 
These states are frequently referred to as the Page mode (left) and the Sievers- 
Takeno mode (right). Adapted from [5]. 




Figure 3: Displacements of DBs on a two-dimensional lattice (j9|) with k = 0.05, 
all velocities equal to zero. (A) ft b = 1.188; (B) fl b = 1.207; (C) Q, b = 1.319. 
From [6]. 

1.2 Prom classical to quantum 

A natural question is what remains of discrete breathers if the corresponding 
quantum problem is considered [7J |H1 [H] • The many-body Schrodinger equation 
is linear and translationally invariant, therefore all eigenstates must obey the 
Bloch theorem. Thus we cannot expect eigenstates of the Hamiltonian to be 
spatially localized (on the lattice). What is the correspondence between the 
quantum eigenvalue problem and the classical dynamical evolution? 

The concept of tunneling is a possible answer to this puzzle. Naively speak- 
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ing we quantize the family of periodic orbits associated with a discrete breather 
located somewhere on the lattice. Notice that there are as many such fami- 
lies as there are lattice sites. The quantization (e.g., Bohr-Sommerfeld) yields 
some eigenvalues. Since we can perform the same procedure with any family 
of discrete breather periodic orbits which differ only in their location on the 
lattice, we obtain ./V-fold degeneracy for every thus obtained eigenvalue, where 
N stands for the number of lattice sites. Unless we consider the trivial case 
of uncoupled lattice sites, these degeneracies will be lifted. Consequently, we 
will instead obtain bands of states with finite band width. These bands will be 
called quantum breather bands. The inverse tunneling time of a semiclassical 
breather from one site to a neighboring one is a measure of the bandwidth. 

We can then formulate the following expectation: if a classical nonlinear 
Hamiltonian lattice possesses discrete breathers, its quantum counterpart should 
show up with nearly degenerate bands of eigenstates, if the classical limit is con- 
sidered. The number of states in such a band is N, and the eigenfunctions are 
given by Bloch-like superpositions of the semiclassical eigenfunctions obtained 
using the mentioned Bohr-Sommerfeld quantization of the classical periodic or- 
bits. By nearly degenerate we mean that the bandwidth of a quantum breather 
band is much smaller than the spacing between different breather bands and the 
average level spacing in the given energy domain, and the classical limit implies 
large eigenvalues. 

Another property of a quantum breather state is that such a state shows 
up with exponential localization in appropriate correlation functions 10 . This 
approach selects all many-particle bound states, no matter how deep one is in 
the quantum regime. In this sense quantum breather states belong to the class 
of many-particle bound states. 

For large energies and N the density of states becomes large too. What will 
happen to the expected quantum breather bands then? Will the hybridization 
with other non-breather states destroy the particle-like nature of the quantum 
breather, or not? What is the impact of the nonintegr ability of most systems 
allowing for classical breather solutions? Since the quantum case corresponds 
to a quantization of the classical phase space, we could expect that chaotic tra- 
jectories lying nearby classical breather solutions might affect the corresponding 
quantum eigenstates. 

From a computational point of view we are very much restricted in our 
abilities to study quantum breathers. Ideally we would like to study quantum 
properties of a lattice problem in the high energy domain (to make contact with 
classical states) and for large lattices. This is typically impossible, since solving 
the quantum problem amounts to diagonalizing the Hamiltonian matrix with 
rank b N where b is the number of states per site, which should be large to make 
contact with classical dynamics. Thus typically quantum breather states have 
been so far obtained numerically for small one-dimensional systems jTOJ [TTJ [TJ] . 
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Figure 4: Spectrum of the quantum DNLS with 6 = 2 and N — 101. The energy 
eigenvalues are plotted versus the wavenumber of the eigenstate. Adapted from 



2 Quantum breather models 
2.1 The Bose-Hubbard chain 

Let us discuss quantum breathers within the widely used quantum discrete 
nonlinear Schrodinger model (also called Bose-Hubbard model) with the Hamil- 
tonian [13] 



N r 



H — ^ 

i=i 



\-a\a\aiai + C(a\ai + \ + h.c.) 



(10) 



2 

and the commutation relations 

aialn - a f m ai = Si m (11) 

with Si m being the standard Kronecker symbol. This Hamiltonian conserves the 
total number of particles 

For b particles and N sites the number of basis states is 

(6 + JV-l)! 



b\(N-l)\ 



(13) 



For b = there is just one trivial state of an empty lattice. For 6=1 there are 
./V states which correspond to one-boson excitations. These states are similar to 
classical extended wave states. For 6 = 2 the problem is still exactly solvable, 
because it corresponds to a two-body problem on a lattice. A corresponding 
numerical solution is sketched in Fig. |4j Note the wide two-particle continuum, 
and a single band located below. This single band corresponds to quasiparticle 
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states characterized by one single quantum number (related to the wavenumbcr 
q). These states are two-particle bound states. The dispersion of this band is 
given |T3] by 

£' = -Jl + 16C 2 cos 2 (|) . (14) 

Any eigenstate from this two-particle bound state band is characterized by ex- 
ponential localization of correlations, i.e. when represented in some set of basis 
states, the amplitude or overlap with a basis state where the two particles are 
separated by some number of sites is exponentially decreasing with increasing 
separation distance. Note that a compact bound state is obtained for q = ±ir, 
i.e. for these wave numbers basis states with nonzero separation distance do 
not contribute to the eigenstate at all [T4"] . 

2.2 The dimer 

A series of papers was devoted to the properties of the quantum dimer [151 1161 
IT7l [T8] . This system describes the dynamics of bosons fluctuating between two 
sites. The number of bosons is conserved, and together with the conservation 
of energy the system appears to be integrable. Of course, one cannot consider 
spatial localization in such a model. However, a reduced form of the discrete 
translational symmetry - namely the permutational symmetry of the two sites - 
can be imposed. Together with the addition of nonlinear terms in the classical 
equations of motion the dimer allows for classical trajectories which are not 
invariant under permutation. The phase space can be completely analyzed, all 
isolated periodic orbits (IPO) can be found. There appears exactly one bifurca- 
tion on one family of isolated periodic orbits, which leads to the appearance of 
a separatrix in phase space. The separatrix separates three regions - one invari- 
ant and two non-invariant under permutations. The subsequent analysis of the 
quantum dimer demonstrated the existence of pairs of eigenstates with nearly 
equal eigenenergies [15] , The separatrix and the bifurcation in the classical 
phase space can be traced in the spectrum of the quantum dimer |17j . 
The classical Hamiltonian may be written as 

F = *t* 1 + ^* 2 + i((* 1 *i) 2 + (*;* 2 ) 2 )+C(1' 1 * 2 + *;*i) . (15) 

with the equations of motion = idH/d^\ 2 - T ne model conserves the norm 
(or number of particles) B — | v I'i| 2 + |*2| 2 - 

Let us parameterize the phase space of (|15l) with = Ai 2& ,a > A 12 > 0. 
It follows that Ax 2 is time independent and </>i — <p2 + A with A = 0, tt and 
0i,2 = w being also time independent. Solving the algebraic equations for the 
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amplitudes of the IPOs we obtain 



I:Al t2 = X -B, A = 0, u=l+C+^B , (16) 
11:^2 = ^5, A = tt, u = 1-C + ±B , (17) 



III :Al = -B\l± y/l-4C 2 /B 2 j , A = , w = 1 + S . (18) 

IPO III corresponds to two elliptic solutions which break the permutational 
symmetry. IPO III exists for B > B\, with B b — 2C and occurs through 
a bifurcation from IPO I. The corresponding separatrix manifold is uniquely 
defined by the energy of IPO I at a given value of B > B b . This manifold 
separates three regions in phase space - two with symmetry broken solutions, 
each one containing one of the IPOs III, and one with symmetry conserving 
solutions containing the elliptic IPO II. The separatrix manifold itself contains 
the hyperbolic IPO I. For B < B b only two IPOs exist - IPO I and II, with both 
of them being of elliptic character. Remarkably there exist no other IPOs, and 
the mentioned bifurcation and separatrix manifolds are the only ones present in 
the classical phase space of (|l~5j) . 

To conclude the analysis of the classical part, we list the energy properties 
of the different phase space parts separated by the separatrix manifold. First 
it is straightforward to show that the IPOs (|16l) - (fT8|) correspond to maxima, 
minima or saddle points of the energy in the allowed energy interval for a given 
value of B, with no other extrema or saddle points present. It follows 

E 1 = if (IPO I) = B + is 2 + CB , (19) 

E 2 = if (IPO II) = B + is 2 - CB , (20) 

E 3 = if (IPO III) = B + ^B 2 + C 2 . (21) 

For B < B b we have E\ > E 2 (IPO I - maximum, IPO II - minimum). For 
B > B b it follows E 3 > E x > E 2 (IPO III - maxima, IPO I - saddle, IPO 
II - minimum). If B < B b , then all trajectories are symmetry conserving. If 
B > Bb, then trajectories with energies E\ < E < E$ are symmetry breaking, 
and trajectories with E 2 < E < E\ are symmetry conserving. 

The quantum eigenvalue problem amounts to replacing the complex func- 
tions ^, in (|15p by the boson annihilation and creation operators a,a^ with 
the standard commutation relations (to enforce the invariance under the ex- 
change W the substitution has to be done on rewriting *vj/* = 1/2('J"I'* + 



H- 5 + 3 



(a\ai + a\a,2j + - (^(a\ai) 2 + (a 2 a 2 ) 2 ^ + C (a\a 2 + a\a^j . (22) 



Note that % = 1 here, and the eigenvalues b of B = a\a\ + a\a 2 are integers. 
Since B commutes with if we can diagonalize the Hamiltonian in the basis of 
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Figure 5: Eigenvalues versus ordered state number h for symmetric and anti- 
symmetric states (0 < h < 6/2 for both types of states). Parameters: b = 600 
and C = 50. Inset: Density of states versus energy. Adapted from |17) . 

eigenfunctions of B. Each value of 6 corresponds to a subspace of the dimension 
(6 + 1) in the space of eigenfunctions of B. These eigenfunctions are products 
of the number states \n) of each degree of freedom and can be characterized by 
a symbol |n, m) with n bosographicsns in the site 1 and m bosons in the site 2. 
For a given value of b it follows m = b — n. So we can actually label each state 
by just one number n: \n, (6 — n)) = \n). Consequently the eigenvalue problem 
at fixed b amounts to diagonalizing the matrix 



where n,m — 0,1,2,..., b. Notice that the matrix H nm is a symmetric band 
matrix. The additional symmetry H nm = H( b _ n - ) ( b _ m - j is a consequence of the 
permutational symmetry of H. For C = the matrix H nm is diagonal, with 
the property that each eigenvalue is doubly degenerate (except for the state 
\b/2) for even values of b). The classical phase space contains only symmetry 
broken trajectories, with the exception of IPO II and the separatrix with IPO 
I (in fact in this limit the separatrix manifold is nothing but a resonant torus 
containing both IPOs I and II). With the exception of the separatrix manifold, 
all tori break permutational symmetry and come in two groups separated by the 
separatrix. Then quantizing each group will lead to pairs of degenerate eigen- 
values - one from each group. There is a clear correspondence to the spectrum 
of the diagonal (C = 0) matrix H nm . The eigenvalues H 00 = correspond to 
the quantized IPOs III. With increasing n the eigenvalues H nn = H( h _ n ) ( b _ n -) 
correspond to quantized tori further away from the IPO III. Finally the states 
with n = 6/2 for even 6 or n = (6 — l)/2 for odd 6 are tori most close to the sep- 
aratrix. Switching the side diagonals on by increasing C will lead to a splitting 
of all pairs of eigenvalues. In the case of small values of 6 these splittings have 




n = m 



m + 1 
m — 1 
else 



(23) 
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Figure 6: Eigenvalue splittings versus n for b — 150 and C = 10. Solid line - 
numerical result, dashed line - perturbation theory. Inset: Same for b = 600 
and C = 50. Only numerical results are shown. Adapted from [17]. 



no correspondence to classical system properties. However, in the limit of large 
b we enter the semiclassical regime, and due to the integr ability of the system, 
eigenfunctions should correspond to tori in the classical phase space which sat- 
isfy the Einstein-Brillouin-Keller quantization rules. Increasing C from zero will 
lead to a splitting AE n of the eigenvalue doublets of C = 0. In other words, we 
find pairs of eigenvalues, which are related to each other through the symmetry 
of their eigenvectors and (for small enough C) through the small value of the 
splitting. These splittings have been calculated numerically and using perturba- 
tion theory [T51 [T7] . In the limit of large b the splittings are exponentially small 
for energies above the classical separatrix energy (i.e. for classical trajectories 
which are not invariant under permutation). If the eigenenergies are lowered 
below the classical separatrix energy, the splittings grow rapidly up to the mean 
level spacing. 

In Fig. [5] the results of a diagonalization of a system with 600 particles 
(b = 600) is shown [T7]. The inset shows the density of states versus energy, 
which nicely confirms the predicted singularity at the energy of the separatrix 
of the classical counterpart. In order to compute the exponentially small split- 
tings, we may use e.g. a Mathematica routine which allows to choose arbitrary 
values for the precision of computations. Here we chose precision 512. In Fig. [6] 
the numerically computed splittings are compared to perturbation theory re- 
sults. As expected, the splittings become extremely small above the separatrix. 
Consequently these states will follow for long times the dynamics of a classical 
broken symmetry state. 

2.3 The trimer 

The integrability of the dimcr does not allow a study of the influence of chaos 
(i.e. nonintegrability) on the tunneling properties of the mentioned pairs of 
eigenstates. A natural extension of the dimer to a trimer adds a third degree of 
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freedom without adding a new integral of motion. Consequently the trimer is 
nonintegrable. A still comparatively simple numerical quantization of the trimer 
allows to study the behavior of many tunneling states in the large-energy domain 
of the eigenvalue spectrum [115] . 

Similarly to the dimer, the quantum trimer Hamiltonian is represented in 
the form 



Again B = a\a\ + a 2 ci2 + ^3^3 commutes with the Hamiltonian, thus we can 
diagonalize (l24l) in the basis of eigenfunctions of B. For any finite eigenvalue b 
of B the number of states is finite, namely (b + l)(b + 2)/2. Thus the infinite 
dimensional Hilbert space separates into an infinite set of finite dimensional 
subspaces, each subspace containing only vectors with a given eigenvalue b. 
These eigenfunctions are products of the number states |n) of each degree of 
freedom and can be characterized by a symbol |n, m, 1} where we have n bosons 
on site 1, m bosons on site 2, and I bosons on site 3. For a given value b it follows 
that I = b — m — n. So we can actually label each state by just two numbers 
(n,m): \n,m, (b — n — m)) = \n,m). Note that the third site added to the dimer 
is different from the first two sites. There is no boson-boson interaction on this 
site. Thus site 3 serves simply as a boson reservoir for the dimer. Dimer bosons 
may now fluctuate from the dimer to the reservoir. The trimer has the same 
permutational symmetry as the dimer. 

The matrix elements of (|24|) between states from different b subspaces van- 
ish. Thus for any given b the task amounts to diagonalizing a finite dimensional 
matrix. The matrix has a tridiagonal block structure, with each diagonal block 
being a dimer matrix (|23[) . The nonzero off-diagonal blocks contain interaction 
terms proportional to 6. We consider symmetric \^) s and antisymmetric \^} a 
states. The structure of the corresponding symmetric and antisymmetric de- 
compositions of H is similar to H itself. In the following we will present results 
for b = 40. We will also drop the first two terms of the RHS in (|2"4")l , because 
these only lead to a shift of the energy spectrum. Since we evaluate the matrix 
elements explicitly, we need only a few seconds to obtain all eigenvalues and 
eigenvectors with the help of standard Fortran routines. In Fig. [7] we plot a 
part of the energy spectrum as a function of 5 for C = 2 [19]. As discussed 
above, the Hamiltonian decomposes into noninteracting blocks for 5 — 0, each 
block corresponding to a dimer with a boson number between and b. For 5 =/= 
the block-block interaction leads to typical features in the spectrum, like, e.g., 
avoided crossings. The full quantum energy spectrum extends roughly over 10 3 , 
leading to an averaged spacing of order 10°. Also the upper third of the spec- 
trum is diluted compared to the lower two thirds. The correspondence to the 
classical model is obtained with the use of the transformation E c i = E qm /b 2 + 1 
and for parameters C/b and S/b (the classical value for B is B = 1). 

The main result of this computation so far is that tunneling pairs of eigen- 
states of the dimer persist in the nonintegrable regime 5^0. However at certain 
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Y + 2( a i fll + a 2 a2 + a 3 a3 ) + 2 ( a i ai ) 2 + (4«2) 2 
+C(a\a 2 + a\ai) + S(a\a 3 + a\ai + a^a 3 + a^a 2 ) 



(24) 
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Figure 7: A part of the eigenenergy spectrum of the quantum trimer as a func- 
tion of 6 with b = 40 and C = 2. Lines connect data points for a given state. 
Solid lines - symmetric eigenstates; thick dashed lines - antisymmetric eigen- 
states. Adapted from [19] . 
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Figure 8: Level splitting versus S for a level pair as described in the text. Solid 
line - numerical result. Dashed line - semiclassical approximation. Filled circles 
- location of wave function analysis in Fig. [9] Adapted from |20) . 

pair-dependent values of 5 a pair breaks up. From the plot in Fig. [7] we cannot 
judge how the pair splittings behave. In Fig[8]we plot the pair splitting of the 
pair which has energy ss 342 at 5 = [20]. Denote with x,y,z the eigenvalues 
of the site number operators 711,712,713. We may consider the quantum states 
of the trimer at S = when z is a good quantum number and then follow the 
evolution of these states with increasing S. The state for S = can be traced 
back to C — and be thus characterized in addition by x and y. The chosen 
pair states are then characterized by x — 26(0), y — 0(26) and z — 14 for 
C = S = 0. Note that this pair survives approximately 30 avoided crossings 
before it is finally destroyed at coupling strength S « 2.67 as seen in Fig. [71 
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Figure 9: Contour plot of the logarithm of the symmetric eigenstate of 
the chosen tunneling pair (cf. Fig. [7]) for five different values of 6 = 
0, 0.3, 0.636, 1.0, 1.8 (their location is indicated by filled circles in Fig. [8]). (a): 
three equidistant grid lines are used; (b-e): ten grid lines are used. Minimum 
value of squared wave function is 10 -30 , maximum value is about 1. Adapted 
from |20]. 

From Fig. [8] we find that the splitting rapidly increases gaining about eight 
orders of magnitude when 5 changes from to slightly above 0.5. Then this 
rapid but nevertheless smooth rise is interrupted by very sharp spikes when 
the splitting AE rises by several orders of magnitude with S changing by mere 
percents and then abruptly changes in the opposite direction sometimes even 
overshooting its pre-spike value. Such spikes, some larger, some smaller, repeat 
with increasing S until the splitting value approaches the mean level spacing of 
order one. Only then one may say that the pair is destroyed since it can be 
hardly distinguished among the other trimer levels. 

Another observation is presented in Fig. EH _2D] . We plot the intensity distri- 
bution of the logarithm of the squared symmetric wave function of our chosen 
pair for five different values of 6 = , 0.3 , 0.636 , 1.0 , 1.8 (their locations are 
indicated by filled circles in Fig. HJ). We use the eigenstates of B as basis states. 
They can be represented as \x, y, z > where x, y, z are the particle numbers on 
sites 1, 2, 3, respectively. Due to the commutation of B with H two site occu- 
pation numbers are enough if the total particle number is fixed. Thus the final 
encoding of states (for a given value of b) can be chosen as \x, z). The abscissa 
in Fig. E2 is x and the ordinate is z. Thus the intensity plots provide us with 
information about the order of particle flow in the course the tunneling process. 
For 6 = (Fig. EH a)) the only possibility for the 26 particles on site 1 is to di- 
rectly tunnel to site 2. Site 3 is decoupled with its 14 particles not participating 
in the process. The squared wave function takes the form of a compact rim 
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Figure 10: Order of tunneling in the trimcr. Filled large circles - sites 1 and 
2, filled small circle - site 3. Arrows indicate direction of transfer of particles. 
Adapted from [20] . 



in the (x, z) plane which is parallel to the x axis. Nonzero values of the wave 
function are observed only on the rim. This direct tunneling has been described 
in chapter \2. 2 1 When switching on some nonzero coupling to the third site, the 
particle number on the dimer (sites 1,2) is not conserved anymore. The third 
site serves as a particle reservoir which is able either to collect particles from 
or supply particles to the dimer. This coupling will allow for nonzero values of 
the wave function away from the rim. But most importantly, it will change the 
shape of the rim. We observe that the rim is bended down to smaller z values 
with increasing 5. That implies that the order of tunneling (when, e.g., going 
from large to small x values) is as follows: first, some particles tunnel from site 
1 to site 2 and simultaneously from site 3 to site 2 (Tig.fTOla)). Afterwards par- 
ticles flow from site 1 to both sites 2 and 3 (Fig. [TU^b)). With increasing 5 the 
structure of the wave function intensity becomes more and more complex, pos- 
sibly revealing information about the classical phase space flow structure. Thus 
we observe three intriguing features. First, the tunneling splitting increases by 
eight orders of magnitude when S increases from zero to 0.5. This seems to be 
unexpected, since at those values perturbation theory in 5 should be applicable 
(at least Fig. 7 indicates that this should be true for the levels themselves). The 
semiclassical explanation of this result was obtained in |20j . 
The second observation is that the tunneling begins with a flow of particles 
from the bath (site 3) directly to the empty site which is to be filled (with si- 
multaneous flow from the filled dimer site to the empty one) . At the end of the 
tunneling process the initially filled dimer site is giving particles back to the 
bath site. Again this is an unexpected result, since it implies that the particle 
number on the dimer is increasing during the tunneling, which seems to decrease 
the tunneling probability, according to the results for an isolated dimer. These 
first two results are closely connected (see [20] for a detailed explanation). 
The third result concerns the resonant structure on top of the smooth variation 
in Fig. [8] The resonant enhancements and suppressions of tunneling are related 
to avoided crossings. Their presence implies that a fine tuning of the system 
parameters may strongly suppress or enhance tunneling which may be useful 
for spectroscopic devices. In Fig. [TT]we show the four various possibilities of 
avoided crossings between a pair and a single level and between two pairs, and 
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Figure 11: Level splitting variation at avoided crossings. Inset: Variation of 
individual eigenvalues participating in the avoided crossing. Solid lines - sym- 
metric eigenstates, dashed lines - antisymmetric eigenstates. Adapted from [20] . 



the schematic outcome for the tunneling splitting [2U] . If the interaction to fur- 
ther more distant states in the spectrum is added, the tunneling splitting can 
become exactly zero |21j for some specific value of the control parameter. In 
such a rare situation the tunneling is suppressed for all times. 

2.4 Quantum roto-breathers 

When discussing classical breather solutions we have been touching some aspects 
of roto-breathers, including their property of being not invariant under time 
reversal symmetry. In a recent study Dorignac et al have performed [22] an 
analysis of the corresponding quantum roto-breather properties in a dimer with 
the Hamiltonian 

H = {— + a(l - cosxj)! + e(l - cos(xi - x 2 )) . (25) 
»=i ^ J 

The classical roto-breather solution consists of one pendulum rotating and the 
other oscillating with a given period Tj,. Since the model has two symme- 
tries - permutation of the indices and time-reversal symmetry - which may be 
both broken by classical trajectories, the irreducible representations of quantum 
eigenstates contain four symmetry sectors (with possible combinations of sym- 
metric or antisymmetric states with respect to the two symmetry operations). 
Consequently, a quantum roto-breather state is belonging to a quadruplet of 
weakly split states rather than to a pair as discussed above. The schematic 
representation of the appearance of such a quadruplet is shown in Fig. [12] [22 . 
The obtained quadruplet has an additional fine structure as compared to the 
tunneling pair of the above considered dimer and trimer. The four levels in 
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Figure 12: Schematic representation of the sum of two pendula spectra. Straight 
solid arrows indicate the levels to be added and dashed arrows the symmetric 
(permutation) operation. The result is indicated in the global spectrum by a 
curved arrow. The construction of the quantum roto-breather state is explicitly 
represented. Adapted from [2"2"j . 




Figure 13: Dependence of different splittings of a quadruplet on e. Only three of 
them have been displayed, each being associated with a given tunneling process. 
Adapted from [22] . 



the quadruplet define three characteristic tunneling processes. Two of them are 
energy or momentum transfer from one pendulum to the other one, while the 
third one corresponds to total momentum reversal (which restores time reversal 
symmetry). The dependence of the corresponding tunneling rates on the cou- 
pling e is shown for a specific quadruplet from |22] in Fig. 1131 For very weak 
coupling e <C 1 the fastest tunneling process will be momentum reversal, since 
tunneling between the pendula is blocked. However as soon as the coupling is 
increased, the momentum reversal turns into the slowest process, with breather 
tunneling from one pendulum to the other one being orders of magnitude faster. 
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Note that again resonant features on these splitting curves are observed, which 
are related to avoided crossings. 



2.5 Large lattices with fluctuating numbers of quanta 

A number of publications are devoted to the properties of quantum breathers 
in chains and two-dimensional lattices of coupled anharmonic oscillators. For 
the respective one-dimensional case, the Hamiltonian is given by 



~pl + V{x n ) + W(x n - X n -i) 



(26) 



Here V(x) = \x 2 + \v±x^ (or similar) and the nearest neighbour coupling 
W(x) — \Cx 2 . The classical version of such models conserves only the energy, 
but not any equivalent of a norm. Therefore, no matter whether one uses 
creation and annihilation operators of the harmonic oscillator [23 , or similar 
operators which diagonalize the single anharmonic oscillator problem |24j . the 
resulting Hamiltonian matrix will not commute with the corresponding number 
operator. Calculations will typically be restricted to 4-6 quanta, and lattice 
sizes of the order of 30 for d = 1, 13 x 13 for d = 2 [23]. With these parameters 
one can calculate properties of quantum breather states, which correspond to 
typically two quanta which are bound (with unavoidable states with different 
number of quanta, contributing as well). For large enough a complete gap 
opens between the two-quanta continuum and quantum breather states [101 123j 
(FigsJ14l and I15p ■ When decreasing the anharmonic constant V4, Proville found, 
that the gap closes for certain wave numbers, but persists for others, becoming 
a pseudogap [25] [23] (FigOU). 

Involved calculations of the dynamical structure factor (e.g. available by 
neutron scattering in crystals) have shown, that signatures of quantum breathers 
are imprinted in these integral characteristics of the underlying lattice dynamics 
p~0| [24] , yet the working out of these differences may become a subtle task (see 
Fig [H] for an example). 

Finally, Fleurov et al [26J estimated the influence of the tunnel splitting of 
a dimer, when embedded in an infinite chain. This situation is close to the 
tunneling of a very localized DB, so that the nonlinearity (interaction between 
bosons) can be taken into account only on the two sites, which participate in 
the tunneling, while the nonlinearity can be neglected on all other sites. Using 
path integral techniques, the computed tunneling splitting has been shown to 
become smaller as compared to the case of an isolated dimer. This is due to 
the fact, that a DB in an infinite chain has a core and a localized tail. That tail 
has to be carried through the tunneling process as well, and in analogy with a 
single particle tunneling in a double well, the tail increases the effective mass of 
such a particle. Consequently the exponential tail of a DB in an infinite chain 
tends to decrease its ability to perform quantum tunneling motion, yet it never 
leads to a full suppression of tunneling [26] . 
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Figure 14: Eigenspectrum of a chain with 33 sites for parameters (a) C = 0.05, 
«4 = 0.2, and (b) C = 0.05, va = 0.02. The inserts show magnifications of the 
fundamental branch (left) and overtone region (right). The quantum breather 
branch is marked by (2), and the two-phonon band by (11). Adapted from [23]. 



3 Quantum breather properties 
3.1 Evolution of quantum localized states 

Suppose that we initially excite only one site in the trimer from above. If this 
initial state has strong enough overlap with tunneling pair eigenstates, its evo- 
lution in time should show distinct properties as compared to the case when the 
overlap is vanishing, or when there are simply no tunneling pair states available. 
Several results have been reported. First, a quantum echo was observed in |19) 
by calculating the survival probability of the initial state as a function of time. 
That quantity measures the probability to find the system in the initial state 
at later times. If the initial state has strong overlap with many eigenstates, it 
is expected to quickly decohere into these different states. Yet, if a substan- 
tial overlap with quantum breathers takes place, the survival probability first 
rapidly decays to zero, but echoes up after regular time intervals (FigfTTl left 
plot). If one simply measures the dependence of the number of quanta, then a 
similar situation will show up with a very slow beating of the occupation num- 
bers in time, if the overlap of the initial state and a tunneling pair is strong 
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Figure 15: Eigenspectrum of a lattice with 13 x 13 sites for parameters (top) 
C = 0.025, v A = 0.1, and (bottom) C = 0.025, v 4 = 0.025. Left plots - spectra 
over the whole Brilloin zone. Right plots - profiles of the spectra along the 
direction [11]. The insets show the magnifications of the phonon branch (left) 
and the quantum breather energy region (right). Adapted from |23j . 




Figure 16: A 3D plot of the inelastic structure factor S(q, w) as a function of the 
dimcnsionless energy transfer < lu < 3 and the scalar product of the transfer 
momentum q and the polarization u. Adapted from |24j . 



[2TJ[27] (see FigOU right plot). 

Suppose we have a large lattice, and put initially many quanta on one site. 
Then any tunneling of this packet as a whole will occur on very long time scales. 
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Figure 17: Left plot: Survival probability of the initial state |\&o) = 1 20 + 
v, 0, 20— v). v = (a) -6, (b) 0, and (c) 6. Insets: spectral intensity of the initial 
state | 'Jo)- Filled circles - symmetric eigenstates; open circles - antisymmetric 
eigenstates. Right plot: Time evolution of expectation values of the number of 
bosons at each site of the trimer for different initial states l^o) = 1 20 + ^ , 0, 20 — 
v). v = (a) -6, (b) 0, and (c) 6. Adapted from gT]- 



On time scales much shorter, we may describe the excitation as a classical 
discrete breather state plus a small perturbation. Treating that perturbation 
quantum mechanically, one could expect that the time-periodic DB acts as a 
constant source of quantum radiation for the quantized phonon field. It turns 
out to be impossible, for very much the same reasons as in the purely classical 
treatment (see [30 ). This result implies, that there is almost no other source of 
decay for a localized initial state in a quantum lattice, but to slowly tunnel as a 
whole along the lattice, if nonlinearities allow for the formation of exact classical 
DB states [31]. Numerical calculations for such a case, but with few quanta, were 
performed by Proville [32], and, similar to the above trimer discussion, showed 
that if quantum breather states exist in the system, then localized excitations 
stay localized for times which are much longer than the typical phonon diffusion 
times in the absence of anharmonicity. 
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Figure 18: (a) Energy splitting and (b) correlation function vs. energy of the 
eigenstates of the dimer (open circles, symmetric eigenstates; solid circles, an- 
tisymmetric eigenstates) . The vertical dashed line marks the energy threshold 
for appearance of QB states. The thin solid line in (a) is a guide for the eye, 
whereas in (b) it is the estimation using eq. (J30j) . Here b — 40 and C = 2. 



3.2 Splitting and correlations 

QBs are nearly degenerate eigenstates. For the dimer and the trimer, they come 
in symmetric-antisymmetric pairs. So one may compute the nearest neighbor 
energy spacing (tunneling splitting) between pairs of symmetric- antisymmetric 
eigenstates in order to identify QBs. Since QBs correspond to classical orbits 
that are characterized by energy localization, they may be identified by defin- 
ing correlation functions. For large lattices it has been shown that QBs have 
exponentially localized correlation functions, in full analogy to their classical 
counterparts. 

For the dimer and the trimer, the correlation functions may be defined as 
follows: 

2) = (n^)^ 1) = (h\)^ (27) 

where n, = afa,, and (A)^ — (Xn\A\x^)i {\Xn)} being the set of eigenstates of 
the system. The ratio < / M (l, 2)// M (l, 1) < 1 measures the site correlation of 
quanta: it is small when quanta are site-correlated (i.e. when many quanta are 
located on one site there are almost none on the other one) and close to unity 
otherwise. 

For the dimer case, the relation b — m + ni , leads to 

f tl {l,2) = b{h 1 ) ll -{hl) tl . (28) 

In Fig. [TBI- left we show the energy splitting and the correlation function of the 
eigenstates. We see that beyond a threshold (dashed line), the splitting drops 
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exponentially fast with energy. The corresponding pairs of eigenstates, which 
are tunneling pairs, are site correlated. Thus they are QBs. Their correlation 
functions show a fast decrease for energies above the threshold. In these states 
many quanta are localized on one site of the dimer and the tunneling time 
of such an excitation from one site to the other (given by the inverse energy 
splitting between the eigenstates of the pair) is exponentially large. As shown 
in Ref. [17], this energy threshold is close to the threshold for the existence of 
DBs in the corresponding classical model. 

QBs are close to symmetric (S) and antisymmetric (A) eigenstates of the 
C = case given by 

\ni,n 2 ) s ,A = -^=(|ni,n 2 ) ± |ra 2 ,rai)), (29) 

with n\ t 2 3> ^2,1- So we may estimate the dependence of the correlation func- 
tions on n± using the eigenstates f|29[) and b = n\ + n 2 - The result is: 

/OA = 2m(6-m) 

/(l,l) ni nl + (b- ni )^ {6V > 

where we note that it is equal to unity when n\ = 6/2, and vanishes when 
ni = 0, b. 

Using the relation between the eigenenergy e of the C = case and the 
number ri\ (= 1, 2, . . . , b) 

*»x = f + [»? + (& -«l) a ], (31) 

one may obtain the energy dependence of the correlation function (|30p . which 
is plotted in Fig. [T8l-b (thin solid line). We can see that beyond the energy 
threshold for appearance of QBs, the numerical results are close to the estima- 
tion flip. 



3.3 Entanglement 

QBs may also be differentiated from other quantum states when measuring the 
degree of entanglement [33] EH]- For the dimer and the trimer the degree of 
entanglement in the eigenstates may be measured by minimizing the distance 
of a given state to the space of product states of the dimer part (expanded by 
the product basis {|ni) (§5 ^2)}), which depends on the largest eigenvalue of the 
corresponding reduced density matrix [35] ESI [37] [55] : 

N 

A = (x«i-™2 - /«i5n 2 ) 2 , (32) 

where for the case of the dimer Xni,n 3 = ( n ii n 2|x), and for the trimer Xni,n 2 — 
(ni,ri2, (77.3 = b — ri\ — n. 2 )|x)- The functions f ni and g n2 are such that A is 
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Figure 19: Left panel: (a) Entanglement of the eigenstates and (b) energy split- 
ting as a function of energy in the dimer (open circles, symmetric eigenstates; 
solid circles, antisymmetric eigenstates). The vertical dashed line marks the 
energy threshold for appearance of QB states. Here b — 40, and C — 2. Right 
panel: The density of the symmetric eigenstates marked by labeled arrows in 
left panel-a: (a) S-0 (arrow 1), (b) S-7 (arrow 2), (c) S-9 (arrow 3), (d) S-19 
(arrow 4). 



minimum |29) . A measures how far a given eigenstate of the system is from 
being a product of single-site states, and has values < A < 1. This measure 
has a direct relation to the distance of a given eigenstate from a possible one 
obtained after performing a Hartree approximation |35] . 

For the dimer, since QB states are close to eigenstates of the C — case 

| X ) Q B^(M)±|0,n», (33) 

with n <b, one expects that the degree of entanglement in QB states is similar 
to the degree of entanglement in such states. Since only two basis states are 
involved, it can not be a state of maximum entanglement. For C = the 
eigenstates of the system are the basis states given by eq. (|29|) . where for 
m = ri2 it follows that A = 0, and for n\ / ri2 (which includes the state in eq. 
[33]) A = 0.5. 

In previous works in a similar quantum dimer model |33[ I34j , it was shown 
that at the energy threshold for appearance of QB states the entanglement (in 
this case measured in a different way) becomes maximum and then decreases 
with energy. From this, and the above reasoning, we expect that QB states 
show decreasing entanglement A with energy, tending to 0.5. Results in left 
panel-a of Fig. [19] agree with this expectation. 

For C = 0, the entanglement has the values and 0.5 corresponding to the 
basis states 1 6/2, 6/2) and \n,b — n) (n / b/2) with equal and distinct number 
of quanta at each site respectively. When C > 0, the eigenstates become linear 
superpositions of the basis states and the entanglement rises, being larger as 
long as more basis states are involved in building up an eigenstate. This can be 
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seen in the right panel of Fig. [T^l where we plot the density p(ni,n 2 = b — ri\) = 
\{n\, b — ni\x)\ 2 = p{n\) of four symmetric eigenstates marked by labeled arrows 
in the left panel: The low-energy eigenstate marked by the arrow 1 consists 
mainly of one basis state: |6/2,6/2), as seen in the right panel-a, hence the 
entanglement is relatively small. When going up in energy the entanglement in 
the eigenstates quickly increases, becoming maximum at the energy threshold, 
and then decreases. An eigenstate just before the threshold like the one marked 
by the arrow 2 in the left panel involves many basis states fulfilling n\ + n% = 
b = 40 (right panel-b), hence the entanglement is large. However, for a QB 
state lying in the energy region beyond the threshold, like the one marked by 
the arrow 3 in the left panel, the number of involved basis states, and thus 
the entanglement starts to decrease (right panel-c). Finally, in high-energy 
eigenstates like the one marked by the arrow 4 in the left panel, which has the 
form shown in eq. (|33[) (right panel-d), the entanglement is even smaller and 
gets close to 0.5 as expected. 

From the above results we see that by measuring entanglement one may 
gain information not only about the energy threshold for existence of QBs (also 
visible when measuring the energy splitting and correlation function), but also 
about how many basis states overlap strongly with the eigenstate under con- 
sideration. We also computed the von Neumann entropy [38] . which is another 
standard measure of entanglement, and the results were consistent with those 
discussed above. 

4 Quantum edge-localized states 

Most of the studies about QBs in large lattices (with few bosons) were done 
considering a system with periodic boundary conditions, and thus, translational 
invariant. However, ususally real systems have to be modeled with open bound- 
ary conditions. Hence it is natural to wonder what hapens with QBs when the 
lattice has finite size, and therefore, no translational invariance. 

In the classical case, it has been shown that in finite nonlinear lattices, 
the breaking of the translational symmetry may lead to the formation of so 
called nonlinear edge states. These are excitations which are localized at the 
edges of the lattice and they have been studied, in particular, in nonlinear 
optics experiments employing optical waveguides, being coined with the name 
of discrete surface solitons [39] . In particular, it is well known that the discrete 
nonlinear Schrodinger equation (DNLS) has time-periodic solutions localized at 
the edge of the lattice [101 E]- It is therefore expected that the large-boson 
limit of the open-boundary Bose-Hubbard model will show eigenstates in which 
the bosons are localized at the edge. 

Numerical studies by Pouthier were done to answer the edge-localization 
question in a lattice with a few bosons [42] , where the mean- field approximation 
(DNLS equation) cannot a priori be expected to provide the correct intuition. 
The answer turns out to be subtle — this phenomenon is not present for the 
case of two particles, but appears when the particle number is three or more 
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Figure 20: Spatial profile of site occupancies for several eigenstates of the three- 
boson Bose-Hubbard chain. The index v counts the eigenstates from the lowest- 
energy one (y = 1). Inset shows the same plot in semilog scale, the linear 
behavior indicating exponential localization in the edge states. This image was 
taken from Ref. [43] 



[42] . Further studies focused on the energy spectrum and eigenstates (Fig. [20)) 
were done by Pinto et al [33], where degenerate perturbation theory allowed to 
explain why edge-localized states exist only if the number of particles is three 
or more [43] . 

5 Quantum breathers with fermions 

Advances in experimental techniques of manipulation of ultracold atoms in op- 
tical lattices make it feasible to explore the physics of few-body interactions. 
Systems with few quantum particles on lattices have new unexpected features 
as compared to the condensed matter case of many-body interactions, where 
excitation energies are typically small compared to the Fermi energy. Therefore 
it is of interest to study binding properties of fermionic pairs with total spin 
zero, as recently presented in Ref. [44]. We use the extended Hubbard model, 
which contains two interaction scales - the on site interaction U and the near- 
est neighbour intersite interaction V. The nonlocal interaction V is added in 
condensed matter physics to emulate remnants of the Coulomb interaction due 
to non-perfect screening of electronic charges. For fermionic ultracold atoms or 
molecules with magnetic or electric dipole-dipole interactions, it can be tuned 
with respect to the local interaction U by modifying the trap geometry of a con- 
densate, additional external dc electric helds, combinations with fast rotating 
external fields, etc (for a review and relevant references see |45|). 

Consider a one-dimensional lattice with / sites and periodic boundary condi- 
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tions described by the extended Hubbard model with the following Hamiltonian: 

H = H + Hu + H v , (34) 

where 

Ho = ~ ^Ahj-^a + %+l,<r) , (35) 
Hu = -U ^n i)t %,4. , %, CT = a+ a a ha , (36) 

3 

H v = -V^2 njfij+i , h,j = nj, t + n hi . (37) 
i 

Ho describes the nearest-neighbor hopping of fermions along the lattice. Here 
the symbols a =t, 4- stand for a fermion with spin up or down. H\j describes 
the onsite interaction between the particles, and Hy the intersite interaction of 
fermions located at adjacent sites. a+ CT and the fcrmionic creation and 

annihilation operators satisfying the corresponding anticommutation relations: 
{a^ CT ,a/,a'} = Sj^Sa-.cr', {a+^a^,} = {a h<7 , ai,a' } = 0. Note that throughout 
this work we consider U and V positive, which leads to bound states located 
below the two-particle continuum. A change of the sign of U, V will simply swap 
the energies. 

To observe the fermionic character of the considered states, any two-particle 
number state is generated from the vacuum \0) by first creating a particle with 
spin down, and then a particle with spin up: e.g. a^a^JO) creates a particle 
with spin down on site 1 and one with spin up on site 2, while a^a^JO) creates 
both particles with spin down and up on site 2. 

Due to periodic boundary conditions the Hamiltonian (IMl) commutes also 
with the translation operator T, which shifts all lattice indices by one. It has 
eigenvalues r = exp(ik), with Bloch wave number k = 2^ and v = 0, 1, 2, /— 
1. 

For the case of having only one fermion (either spin up or spin down) in 
the lattice (n = 1), a number state has the form \j) = a^ a \0). The interaction 

terms Hjj and Hy do not contribute. For a given wave number fc, the eigenstatc 
to (fM)) is therefore given by: 

i^^EgrV). (38) 

v J s=l 

The corresponding eigenenergy 

e k = -2cos(fc). (39) 

For two particles, the number state method involves N s = f 2 basis states, 
which is the number of ways one can distribute two fermions with opposite 
spins over the / sites including possible double occupancy of a site. Below we 
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Figure 21: (Color online) Energy spectrum of the two fermion states. The 
eigenvalues are plotted as a function of the wave number k. Here U = 1 , V — 
1, / = 101. Symbols are from analytical derivation, lines are the result of 
numerical diagonalization. The arrows indicate the location of the critical wave 
numbers (see text). Adapted from [31] . 



consider only cases of odd / for simplicity. Extension to even values of / is 
straightforward. The details of the calculations can be found in [44] , 

In Fig. [21]we show the energy spectrum of the Hamiltonian matrix obtained 
by numerical diagonalization for the interaction parameters U = 2 and V — 
2 and / = 101. At U = and V — 0, the spectrum is given by the two 
fermion continuum, whose eigenstates are characterized by the two fermions 
independently moving along the lattice. In this case the eigenenergies are the 
sum of the two single-particle energies: 

E° kl M =-2[cos(fci)+cos(fc 2 )] 1 (40) 

with ki t 2 = nvia/if + 1) > v i,n = !)•••»/■ The Bloch wave number k = k\ + &2 
mod 2ir. Therefore, if k — ±tt, the continuum degenerates into points. The con- 
tinuum is bounded by the hull curves h±{k) — ±4 cos |. The same two-particle 
continuum is still observed in Fig. [5T]for nonzero interaction. However, in addi- 
tion to the continuum, we observe one, two or three bound states dropping out 
of the continuum, which depends on the wave number. For any nonzero U and 
V, all three bound states drop out of the continuum at k — ±tt. One of them 
stays bounded for all values of k. The two other ones merge with the continuum 
at some critical value of upon approaching k — as observed in Fig|2"Tl Note 
that for k = ±tt and U = V, all three bound states are degenerate. 

Upon increasing U and V, we observe that a second bound state band sep- 
arates from the continuum for all k fFigl^j). At the same time, when U ^ V, 
the degeneracy at k = ±7r is reduced to two. 

Finally, for even larger values of U and V, all three bound state bands 
completely separate from the continuum (FigJ23p. 
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Figure 22: (Color online) Energy spectrum of the two fermion states. The 
eigenvalues are plotted as a function of the wave number k. Here U = 4 , V — 
3 , / = 101. The symbols are from analytical derivation, lines are the results of 
numerical diagonalization. The arrow indicates the location of the critical wave 
number (see text). Adapted from [4"4"j . 



5.1 Symmetric and antisymmetric state representation 

In order to obtain analytical estimates on the properties of the observed bound 
states, we use the fact that the Hamiltonian for a two fermion state is invariant 
under flipping the spins of both particles. We define symmetric basis states 

|$ J> ) = -^(|$ i , + > + |^,_» (41) 

and antisymmetric states 

|* i)B ) = - !*;,-» • (42) 

Note that |$i) is a symmetric state as well. 

5.2 Antisymmetric bound states 

The antisymmetric states exclude double occupation. Therefore the spectrum 
is identical with the one of two spinless fermions [13] . Following the derivations 
in [13] we find that the antisymmetric bound state, if it exists, has an energy 

E«(k) = -(V+^cos 2 (t)). (43) 

This result is valid as long as the bound state energy stays outside of the contin- 
uum. The critical value of k at which validity is lost, is obtained by requesting 
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Figure 23: (Color online) Energy spectrum for the two fermion states. The 
eigenvalues are plotted as a function of the wave number k. Here U = 8 , V — 
8 , / = 101. The symbols are from analytical derivation, lines are the results of 
numerical diagonalization. Adapted from [44] . 



I^WI = l^±(^)l- K follows V — 2cos(|). Therefore the antisymmetric bound 
state merges with the continuum at a critical wave number 

k a c = 2arccos(-^) , (44) 

setting a critical length scale A" = |j. For V — 1 it follows k^/n « 0.667 (see 
FigHU. 

The equation (|43[) is in excellent agreement with the numerical data in Figs. 
1211221231 fcf. open triangles). We also note, that the antisymmetric bound state 
is located between the two symmetric bound states, which we discuss next. 



5.3 Symmetric bound states 

A bound state can be searched for by assuming an unnormalized eigenvector of 
the form |c, 1, fx, p 2 , p 3 , ...) with \p\ = p < 1. We obtain [44] 

Ef(fc) = -2(p+ i)cosfc/2 . (45) 

The parameter p satisfies a cubic equation 

ap 3 + bp 2 +cp + d = (46) 

with the real coefficients a, b. c and d given by a = 2Fcos(|), b = 4cos 2 (|) — 
UV, c = 2(U + V)cos(^), d = -4cos 2 (|). We plot the results in Figs. I21l22l23l 
(cf. open circles and squares). We obtain excellent agreement. 



30 



At the Brilloin zone edge k — ±tt the cubic equation (f4l)[) is reduced to a 
quadratic one, and can be solved to obtain finally p — > and 



E^(k -> ±tt) 



U , Ef(k -> ±tt) = -V . 



(47) 



In particular we find for k = ±tt that -E| 2 = E%- In addition, if U = V, all 
three bound states degenerate at the zone edge. 

If V = 0, the cubic equation (|4"B]l is reduced to a quadratic one in the whole 
range of k and yields [13] 



Next we determine the critical value of k for which the bound state with 
energy i?f 2 is joining the continuum. Since at this point p = 1, we solve (|46t 
with respect to /c c and find 



setting another critical length scale A* = |f. E.g. for U — V = 1 fc^ 2 /7r ss 0.89, 

in excellent agreement with Fig|2T] For [7 = 4 and = 3 we find k^ 2 /ir « 0.59 
confirming numerical results in Fig |22l 

6 Small Josephson junction networks 

Recent studies of Pinto et al [2H1 ESJ deal with quantum breather excitations in 
two capacitively coupled Josephson junctions. Such systems are currently under 
experimental investigation, being candidates for quantum information process- 
ing, and show remarkably long coherence times up to 100 ns for few quanta 
excitations. The system does not conserve the number of excited quanta, and 
can be best compared with the above Bose-Hubbard trimer. Quantum breather 
signatures are found simultaneously in the spectra (tunneling splittings), corre- 
lation functions, entanglement, and quanta number fluctuations. 

We address the excitation of QBs in a system of two coupled Josephson junc- 
tions in the quantum regime |29j . Josephson junctions are nonlinear devices that 
show macroscopic quantum behavior, and nowadays they can be manipulated 
with high precision, in such a way that the energy flow between coupled junc- 
tions can be resolved in time. 

Josephson junctions behave like anharmonic oscillators, and by lowering the 
temperature one can bring them into the quantum regime, where their quan- 
tization leads to energy levels which are nonequidistant because of the anhar- 
monicity. These levels can be separately excited by using microwaves pulses, 
and the energy distribution between the junctions can be measured in time using 
subsequent pulses |46j . So far these techniques have been used for experiments 
on quantum information processing with Josephson junctions [471 148[ [49], but 
we think that arrays of Josephson junctions in the qunatum regime also might 



£ 2 sl (fc) = -^U 2 + 16cos 2 (fc/2) • 



(48) 




(49) 
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(a) (b) (c) 

Figure 24: (a) Sketch of the two capacitively coupled Josephson junctions, 
(b) Sketch of the washboard potential for a single current-biased Josephson- 
junction. (c) Circuit diagram for two ideal capacitively coupled Josephson junc- 
tions. 



be used as playgrounds for experiments on quantum dynamics of excitations in 
nonlinear lattices. 

The system is sketched in Figl2"4la: two JJs are coupled by a capacitance C c , 
and they are biased by the same current The strength of the coupling due to 
the capacitor is ( — C c /{C C + Cj). The dynamics of a biased Josephson junction 
(JJ) is analogous to the dynamics of a particle with a mass proportional to the 
junction capacitance Cj, moving on a tilted washboard potential 

L%) = -I c -9- cos ip - hp-?-, (50) 

which is sketched in Figl24l-b. Here ip is the phase difference between the macro- 
scopic wave functions in both superconducting electrodes of the junction, If, is 
the bias current, I c is the critical current of the junction, and $o = h/2e the 
flux quantum. When the energy of the particle is large enough to overcome the 
barrier AU (that depends on the bias current lb) it escapes and moves down the 
potential, switching the junction into a resistive state with a nonzero voltage 
proportional to tp across it. Quantization of the system leads to discrete energy 
levels inside the potential wells, which are nonequidistant because of the anhar- 
monicity. Note that even if there is not enough energy to classically overcome 
the barrier, the particle may perform a quantum escape and tunnel outside the 
well, thus switching the junction into the resistive state. Thus each state inside 
the well is characterized by a bias and state-dependent lifetime, or its inverse 
— the escape rate. 

The Hamiltonian of the system is 

H=^ + p- + U( l p 1 ) + U(i P2 ) + ^P 1 P 2 , (51) 
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Figure 25: (a) Energy splitting and (b) correlation function vs. energy of the 
eigenstates of the two-junctions system (open circles, symmetric eigenstates; 
filled circles, antisymmetric eigenstates). The labeled arrows mark the energy 
corresponding to the peak of the spectral intensity in Figl26l-b. d, and f (see 
text). The parameters are 7 = Ib/Ic = 0.945 and £ = 0.1 (22 levels per 
junction). 



where 

m = Cj(l + C)(^\ , (52) 

Pi,2 = ( C ^ + C J){^) (<Pi,2-CM- (53) 

Note that the conjugate momenta Pi 2 are proportional to the charge at the 
nodes of the circuit (which are labeled in FigJU-c). 

In the quantum case the energy eigenvalues and the eigenstates of the system 
were computed and analyzed in 29 . In Fig |25l we show the nearest neighbor 
energy spacing (tunneling splitting) and the correlation function of the eigen- 
states. For this, and all the rest, we used I c = 13.3 [iA, Cj = 4.3 pF, and 
C = 0.1, which are typical values in experiments. We see that in the central 
part of the spectrum the energy splitting becomes small in comparison to the 
average. The corresponding pairs of eigenstates, which are tunneling pairs, are 
site correlated, and thus QBs. In these states many quanta are localized on 
one junction and the tunneling time of such an excitation from one junction 
to the other (given by the inverse energy splitting between the eigenstates of 
the pair) can be exponentially large and depend sensitively on the number of 
quanta excited. 

Note that the tunneling of quanta between the J Js occurs without an obvious 
potential energy barrier being present (the interaction between the junctions is 
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Figure 26: Time evolution of expectation values of the number of quanta at 
each junction (left panels) for different initial excitations with corresponding 
spectral intensities (right panels), (a) and (b): \^q) = |0, 5); (c) and (d): 
l^o) = |0, 19); (e) and (f): |\l>o) = 1 9, 19) . Open circles, symmetric eigenstates; 
filled circles, antisymmetric eigenstates. The energies of the peaks in the spectral 
intensity are marked by labeled arrows in Fig |25l -b (see text). The parameters 
are 7 = Ib/Ic = 0.945 and £ = 0.1 (22 levels per junction). 



only through their momenta). This process has been coined dynamical tun- 
neling [50l |5ll E2] , to distinguish from the usual tunneling through a potential 
barrier. In dynamical tunneling, the barrier — a so-called invariant separatrix 
manifold — is only visible in phase space, where it separates two regions of regu- 
lar classical motion between which the tunneling process takes place. Therefore, 
when referring to the tunneling between the JJs, we implicitly mean that it is 
dynamical. 

The fact that the strongest site correlated eigenstates occur in the central 
part of the energy spectrum may be easily explained as follows: Let N be the 
highest excited state in a single junction, with a corresponding maximum energy 
AU (Figi24]). For two junctions the energy of the system with both junctions 
in the JV-th state is 2AU, which roughly is the width of the full spectrum. 
Thus states of the form \N, 0) and |0, N) that have energy AU are located 
approximately in the middle. 

Having the eigenvalues and eigenstates, we compute the time evolution of 
different initially localized excitations, and the expectation value of the num- 
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n 1 

Figure 27: Contour plot of the logarithm of the density of the asymmetric state 
Ix) = dXif ^) + lxi A ''))/v / 2 as a function of the number of quanta at junctions 1 
and 2 (see text). The parameters are 7 = Ib/Ic = 0.945 and £ = 0.1 (22 levels 
per junction). 



ber of quanta at each junction = (^(i)|fii|\E , (t)). Results are shown in 

Figl2l)R. c, and e. We also compute the spectral intensity 1° = |(x^|*o)| 2 , which 
measures the strength of overlap of the initial state l^o) with the eigenstates. 
Results are shown in Fig |26l b. d, and f, where we can see a peak in each case, 
which corresponds to the arrows in Fig|25l-b. We can see that the initial state 
l^o) = |0, 5) overlaps with eigenstates with an energy splitting between them 
being relatively large and hence the tunneling time of the initially localized 
excitation is short. For the case l^o) = 1 0, 19) QBs are excited: The excita- 
tion overlaps strongly with tunneling pairs of eigenstates in the central part of 
the spectrum, which are site correlated and nearly degenerate. The tunneling 
time of such an excitation is very long, and thus keeps the quanta localized 
on their initial excitation site for corresponding times. Finally the initial state 
l^o) = |9, 19) overlaps with weakly site correlated eigenstates with large energy 
splitting. Hence the tunneling time is short again. 

We computed also the time evolution of the expectation values of the number 
of quanta for initial conditions which are coherent or incoherent (mixtures) 
superpositions of product basis states with equal weights. This is relevant for 
experiments, since it may be hard to excite one junction to a determined state 
but easier to excite several states of the junction at the same time. We used 
coherent superpositions (characterized by well defined states |\l/o)), an d mixtures 
(characterized by their corresponding density operators po), of four basis states 
around the already used initial states: For the state |0,5) we superposed the 
basis states |0,5), |0,6), |0,7), and |0,8), for |0, 19) the basis states |0, 20) , 
|0, 19), |0, 18), and |0,17), and for |9, 19) the basis states |9,20), |9, 19), |9, 18), 
and |9, 17). Both for superposition and mixture of basis states, the results are 
qualitatively similar to those shown in Figj26] Therefore we expect that some 
imprecision in exciting an initial state in the junctions would not affect in a 
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relevant way the results, we may also conclude, that the excitation of QB states 
does not rely on the phase coherence. 

Let us estimate how many quanta should be excited in the junctions in 
order to obtain QBs (tunneling pairs). We compute the density p[n\,fi2) — 

|("-i,"-2|x)| 2 of the asymmetric state \x) = {\x[ S) ) + \xi A) ))/V2, where \xi S ' A) ) 
are the eigenstates belonging to a tunneling pair 29 . In Fig. [57] we show a 
contour plot of the logarithm of the density for the tunneling pair with energy 
marked by the arrow labeled by number two in Fig. [25]- (b). We see that the 
density has its maximum around ri\ = 19 and 712 — which is consistent with 
the result shown in Fig. I2"6l c and d where QBs were excited by using this 
combination of quanta in the junctions. 

7 Experimental realizations 

There is a fast growing amount of experimental and related theoretical work on 
applying the quantum discrete breather concept to many different branches in 
physics, like Bose-Einstein condensates, crystals and molecules, surfaces, and 
others. We will discuss some of these at length, while others will be reviewed 
more briefly. 

7.1 Repulsively bound atom pairs 

Winkler et al [S3] performed experiments with a three-dimensional optical lat- 
tice with initially each site being either not occupied, or being occupied by two 
Rb atoms bound in a pair due to attractive interaction. A magnetic field sweep 
across the Feshbach resonance changes the sign of interaction, turning attraction 
into repulsion. The dynamics of ultracold atoms loaded into the lowest band of 
the optical potential is described by the quantum DNLS model, which is equiv- 
alent to the Bose- Hubbard model (p~0|) . Lifetime measurements have shown, 
that repulsive pairs of Rb atoms have larger lifetimes than pairs of weakly or al- 
most not interacting atoms (Tig|28l) . The two-particle bound states discussed in 
chapter |2. II - the simplest versions of a quantum discrete breather - are the obvi- 
ous explanation of the experimental findings. Indeed, neglecting Landau-Zener 
transitions to higher lying bands in the optical potential, the Bose-Hubbard 
model is justified. The sign of the interaction does not play any role, since it 
only changes quantum discrete breathers from being low-lying to being excited 
states, not affecting their localization properties. The most simple argument of 
why two quanta (or atoms) placed initially close to each other, do not separate 
despite they repel each other, is based on the fact, that if they would do so, 
the (large) interaction energy should be converted into kinetic energy, which is 
restricted to be less than two times the width of the single particle band. In 
other words, repulsively bound atom pairs are a straightforward consequence of 
quantum discrete breather states with two quanta. 

Another sophisticated experimental investigation aimed at measuring the 
quasi-momentum distribution of atom pairs in various regimes by mapping it 
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Figure 28: (a) Repulsive interaction between two atoms sharing a lattice site 
gives rise to an interaction energy U . Breaking up of the pair is suppressed 
owing to the lattice band structure and energy conservation. This is the simplest 
version of a quantum discrete breather, (b) The discrete breather can tunnel 
along the lattice, (c) Long lifetimes of strongly repulsive atoms. The plot shows 
the remaining fraction of pairs for strong interaction (open diamonds) and for 
almost vanishing interaction (filled circles). Adapted from 53 . 



onto a spatial distribution, which was finally measured using standard absorp- 
tion imaging (Fig l29|) . Therefore predictions of such states, which were made 
more than 30 years ago by Ovchinnikov [54], were beautifully confirmed exper- 
imentally with ultracold repulsive atoms. 

7.2 Molecules 

Intramolecular vibrational energy redistribution (IVR) has been a central is- 
sue in the field of chemical physics for many decades. In particular, pathways 
and rates are of importance there, since understanding them allows to describe 
e.g. the dynamics of various chemical reactions, and dissociation processes [55] . 
Spectroscopical studies, where single vibrational quanta are excited, allow to 
measure the frequencies of normal vibrational modes, i.e. to characterize the 
dynamics of a molecule for small amplitude vibrations. These normal modes 
consist of coherent combinations of vibrational excitations of several bonds (or 
rotational groups) in a molecule. However, in order e.g. to dissociate a molecule, 
a many quanta excitation is needed, and nonlinearities will certainly become im- 
portant. It was realized then, that strong vibrational excitations of molecules are 
much better described by so-called local modes, i.e. basically one or few bond 
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Figure 29: (Upper row) Absorption images of the atomic distribution after 
release from the 3D lattice and a subsequent 15 ms time of flight. The horizontal 
and vertical black lines enclose the first Brillouin zone, (upper left) Lattice sites 
are occupied by single atoms; (upper middle) Repulsively bound atom pairs; 
(upper right) Attractively bound atom pairs; (bottom row) Quasi-momentum 
distribution for pairs in one direction as a function of the lattice depth after 
integrating out the other direction, (bottom left) Experiment; (bottom right) 
Numerical calculations. Adapted from |53j . 



vibrational excitations. That transition from normal to local modes remained 
a puzzle for a long time. A practically complete modern theoretical account on 
these issues can be found in a recent monography by Ovchinnikov, Erikhman 
and Pronin 56 . On its most abstract level, the transition from normal to local 
modes is identical with the bifurcation in the dimer model. Thus, local modes 
are essentially discrete breathers or slight perturbations of them. Note, that 
the connection between local modes, breathers and periodic orbits has been 
recently studied by Farantos in the context of large biological molecules [57] . 
Discrete breathers (ILMs) have been theoretically predicted to exist in ionic 
crystals [58] , ways of optical excitation of DBs (ILMs) have been proposed [59] , 
and their possible presence in hydrocarbon structures has been discussed |60) . 

Exciting local modes in molecules with discrete symmetries leads to small 
tunneling splittings of excitation levels [56] , and goes back to the work of Child 
and Lawton |61) , see also a recent comprehensive review by Keshavamurthy [62j 
and references therein. On its most abstract level, this effect is identical with 
the tunneling splitting in the permutationally symmetric dimer model discussed 
in chapter \2. 2 1 

An early example of experimental evidence of discrete breather excitations in 
molecules comes from spectroscopical studies of visible red absorption spectra of 
benzene, naphtalene, and anthracene by Swofford et al [64] ■ The C-H stretching 
vibrations have been excited to the sixth quantum level, and red shifts of the 
lines show, that instead of a delocalized excitation of six bonds to the first 
level (yielding six quanta) , the excitation resides on just one of the six available 
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Figure 30: (a) Absorption spectra of PBLG in chloroform at 293 K (red line, 
helical conformation) and at 260 K (blue line, random coil), (b) Pump-probe 
spectra 600 fs after excita tion under the same conditions. Inset: Decay of 
negative and both positive bands at 293 K. Adapted from [63 . 



bonds. While it can tunnel (as a quantum discrete breather) to the other bonds, 
this tunneling time is a new large time scale in the problem, strongly affecting 
e.g. dissociation rates. 

A recent study of femtosecond infrared pump-probe spectroscopy of the N- 
H mode of a stable a-helix (poly-7-benzyl-L-glutamatc (PBLG)) revealed two 
excited-state absorption bands, which disappear upon unfolding of the helix [63] . 
PBLG forms extremely stable, long a-helices in both helicogenic solvents and 
films grown from these solvents. The monomeric unit of PBLG is a non-natural 
amino-acid with a long side chain that stabilizes the helix. PBLG has served as 
the standard model helix since the very early days of structural investigations 
of proteins. Fig lBUT a) (red line) shows the absorption spectrum of the helix at 
293 K, which is dominated by the strong N-H stretching band at 3290 cm -1 . 
Fig l30f b) (red line) shows the pump-probe response 600 fs after excitation with 
an ultrashort broadband pulse. One negative (3280 cm" 1 ) and two positive 
bands (3160 and 3005 cm -1 ) are observed. If the N-H stretching modes were 
isolated, a negative band associated with bleach and stimulated emission, and a 
positive band associated with excited-state absorption, would be expected. This 



39 




Figure 31: Simulated pump-probe spectrum for 293 K (red line) and 18 K (blue 
line). Inset: schematic of the energy levels. Adapted from [63.. 



is indeed observed here for the unfolded molecule. In contrast, the observation 
of two positive bands for the intact helix rather than just one, is exceptional. 
Edler et al [63] argue, that these features can not be explained due to inten- 
sity dependencies, or Fermi resonances. A consistent explanation is reached 
by assuming that two vibron states are excited, and these vibrons may form 
two different types of bound states, self trapped either on the same site, or on 
neighbouring ones. The latter states originate from the acoustic phonons of the 
helix, which correlate adjacent sites (see also [65]). 

7.3 Crystal surfaces 

Depositing atoms or molecules on crystal surfaces can be controlled experimen- 
tally, and as a result a planar regular two-dimensional lattice structure of the 
deposited material can be obtained. Guyot-Sionnest [66] used Hydrogen to be 
deposited on Si(lll) surfaces. The Si-H bonds can be excited using pum-probe 
techniques with infrared dye lasers. There is substantial interaction between the 
Si-H bonds on the Si(lll) surface. The pump excites one phonon (quantum), 
while the tunable probe frequency finds a substantial red shift of the two-phonon 
excitation, and allows to conclude, that two-phonon bound states are observed. 

Another set of experiments was performed by Jakob [67, 68, 69 . Carbon 
monoxide (CO) was deposited on a Ru(OOl) single crystal surface. The C-0 
stretching modes constitute a two-dimensional array of weakly interacting an- 



40 




1325 1950 1975 2909 2925 2059 3900 3925 3950 3S75 4099 

Wave number cm Wave number cm' 1 



Figure 32: (a) Infrared absorption spectra of the C-0 stretching mode at 30 
K. The corresponding mode of naturally abundant 13 C 16 is displayed in an 
enlarged vertical scale; (b) the overtone band observed at less than twice the 
frequency of the fundamental mode. Adapted from [57] . 
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Figure 33: Frequency shifts of the vibrational bands with temperature: crosses 

- overtone band, dots - fundamental of naturally abundant 13 C 16 0, open circles 

- delocalized fundamental of 12 C 16 0. Adapted from [67] . 



harmonic oscillators with 4.7 A intermolecular distance. Intermolecular coupling 
is provided by means of the electric field of the oscillating dipoles. Experimen- 
tal spectra at 30 K are shown in Fig l32l The one phonon mode frequency is 
at 2031 cm -1 . This has to be compared to the naturally abundant 13 C 16 
frequency at 1941 cm -1 . The corresponding blue shift for the adsorbate is thus 
due to additional stiffness provided by the Ru surface coupling. Excitation of 
two uncorrelated phonons would yield a two phonon continuum at about 4062 
cm . The narrow line observed at 3940 cm -1 can be thus attributed to a 
two-phonon bound state, or a quantum discrete breather excitation. 

The temperature dependence of the line positions also clearly shows, that the 
two-phonon bound state line softens much slower than the line of the one-phonon 
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Figure 34: Structure of the PtCl crystal. One PtCl chain is shown on the left. 
Each Pt atom is coordinated by two ethylenediamine units in a near square 
planar geometry, while CI ions connect pt sites along the chain. The packing 
arrangement of the ID chains and their CIO4 counterions is shown on the right. 
Adapted from [76] . 



delocalized state fFig l33| . This is, among other facts, a strong indication that 
the observed red shift of the overtone line is due to the formation of a localized 
two-phonon bound state, or a (quantum) discrete breather. 

7.3.1 In the bulk of solids 

Vibrational spectra in the overtone or combination region of molecular crystals 
have been studied intensively in the 1970s and 1980s. A pioneering theoretical 
proposal was due to Agranovich, who predicted the existence of two-exciton 
bound states in various molecular crystal materials [70]. Experimental studies 
of infrared absorption spectra for CO2 crystals were conducted by Dows et al 
[7T] and gave evidence of two-phonon bound states. Dressier et al studied the 
slow vibrational relaxation of N2, which also indicates the presence of many- 
phonon bound states [72]. In a remarkable theoretical paper, Bogani calculated 
the spectrum of two phonon excitations in molecular crystals [73], to some 
extent one of the first accurate calculations of quantum discrete breathers. More 
recently Bini et al reconsidered the theory of three-phonon bound states in 
crystal CO2 [74]. While there certainly are many other results worth to be 
mentioned, we recommend reading related chapters in [75], [55]. 

The pioneering studies of Swanson et al |76j have shown that up to seven 
phonons can bind and form a localized state. The system of choice was a PtCl 
based crystal - a halide-bridged mixed-valence transition metal complex, which 
is a model low-dimensional electronic material where the ground states can be 
systematically tuned (with chemistry, doping, pressure, and temperature). It 
is a very strong charge-density wave (CDW) example. The material is a well- 
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Figure 35: Fundamental and overtone spectra of isotopically pure Pt 35 Cl. Mov- 
ing upward in each panel, each x axis is offset by the appropriate integral mul- 
tiple of the 312 cm -1 fundamental frequency. All spectra have been scaled 
vertically to equal peak intensities. Adapted from (76j . 



formed crystal with a homogeneous lattice consisting of quasi-lD chains (see 
Fig |34|) . The CDW ground state consists of alternating Pt +2 and Pt +4 sites 
with a corresponding distortion of the chloride ions towards the Pt +4 site. Res- 
onance Raman spectra were used to probe both ground and photoexcited states. 
They probe the fundamental Cl-Pt-Cl stretch and the progression of many over- 
tones. At low temperatures, the fundamental exhibits a fine structure with up 
to six discrete, well-resolved modes. The analysis of the evolution of the spec- 
tral structure in the overtones was performed for isotopically pure samples, in 
order to avoid exciting localized states due to isotopic disorder. The fundamen- 
tal and overtone spectra for the pure Pt 35 Cl sample are shown in Fig 135] The 
data are presented in a stack plot in which each successive trace is offset along 
the horizontal axis by increasing multiples of the fundamental frequency 312 
cm -1 . Such plots clearly expose the relation of features in the overtone spec- 
trum to multiples of the fundamental peak. The lowest energy dominant feature 
in each trace (marked by vertical lines) demonstrates a strongly increasing an- 
harmonic redshift. Further, at higher overtones, each of these dominant peaks 
recurs, offset by the fundamental frequency, in the next trace above. A simple 
interpretation is that the lowest-energy dominant peaks in the overtone spectra 
correspond to all quanta of vibrational energy localized in approximately one 
PtCl2 unit, while the higher energy peaks correspond to having all quanta but 
one in a localized PtCl2 unit combined with one quanta in the more extended 
fundamental. The schematic process of the energy transfer is shown in Fig |36l 
and has been analyzed theoretically in |77) . 
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Figure 36: A simple picture of a resonant Raman scattering event in the localized 
atomic limit. Large filled circles mark Pt ions, small grey circles mark CI ions. 
Open circles mark the positions of electrons. Adapted from [77] . 



An incoming photon at frequency v is exciting an electron from a Pt +2 to 
a Pt +4 site. The CI ion between them starts oscillating. Finally the electron 
relaxes back to its original position, and releases a photon with frequency v' . 
The energy difference is remaining in a localized vibration. The effect of isotope 
disorder was treated by Kalosakas et al [78]. The experimentally observed red- 
shifts were also theoretically described by Fehske et al [73] using Peierls- Hubbard 
models, and by Wellein et al using the Holstein model [80]. 

Inelastic X-ray and neutron scattering was used by Manley et al [81" ] l82 ] l83] to 
probe phonon dispersion in a-uranium single crystals. Variation of temperature 
showed softening, and the abrupt appearance of a new dynamical mode, without 
a typically observed phase transition. The authors argue that this mode is a 
discrete breather, and forms due to strong electron-phonon interaction. 

Russell and Eilbeck reported evidence for moving breathers in the layered 
crystal muscovite at 300 K [84] . Breathers were created by bombardment of the 
crystal surface with heavy ions. Ejection of atoms at the opposite (shielded) 
crystal surface was attributed to breathers, which were able to carry the vibra- 
tional energy without dispersing over more than 10 unit cells of the crystal. 
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Finally Abrasonis et al [55] reported on anomalous bulk diffusion of inter- 
stitial Nitrogen in steel microcrystals. N ions were deposited in a micron-thick 
layer, and are trapped by local structures, with a characteristic binding energy. 
Ar ion bombardment increases the N mobility at depths far beyond the ion 
penetration depth. The authors see evidence for coherent transfer of vibrational 
energy deep into the bulk of the material. 

8 Conclusions and outlook 

Progress in the understanding of classical discrete breathers evolved for two 
decades, and has significantly improved our understanding of the complexity in 
the dynamics of anharmonic lattice dynamics. The little input it needs to form 
discrete breathers - spatial discreteness and nonlinearity - were demonstrated 
to generate an impressive list of experimental observations in a large variety of 
different physical systems, with length and time scales ranging over many orders 
of magnitude. The quantization of the classical equations of motion opened a 
whole new avenue of research on quantum breathers. For obvious computational 
reasons numerical studies are much more involved in the quantum domain, re- 
stricting our potential of detailed understanding to either small lattices, or few 
participating quanta. Nevertheless, a set of experiments, comparable in size 
to its classical counterpart, has been reviewed where quantum breathers are 
detected. Still we anticipate this only to be the beginning of a much more so- 
phisticated undertaking. We may expect many new qualities to show up in the 
regime of large lattices and many participating quanta. These limits are concep- 
tually closing the gap to quantum interacting many body problems, which are 
at the heart of condensed matter physics, and quantum computing. Yet another 
highly interesting, complicated and debated problem is the combined influence 
of discreteness, nonlinearity, and disorder. Linear disordered wave equations 
allow for Anderson localization. It is a thrilling modern theme to study the 
constructive and destructive interplay of all three ingredients. The best is yet 
to come. 
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